PUPT-1410 
hep-th/9307053 



Light-Cone Quantisation of Matrix Models at c > 1 



Simon Dalley 

Joseph Henry Laboratories 
Princeton University 
Princeton, New Jersey 08544 



ABSTRACT 

The technique of (discretised) hght-cone quantisation, as apphed to matrix 
models of relativistic strings, is reviewed. The case of the c = 2 non-critical 
bosonic string is discussed in some detail to clarify the nature of the continuum 
limit. Further applications for the technique are then outlined. 
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1. Introduction 



Random surface problems appear in many branches of theoretical and mathematical 
physics. A number of them may be reformulated as, or arise from, matrix field theories 
[1, 2, 3, 4, 5, 6], planar diagrams modelling the fluctuating surfaces. The possible physical 
applications, together with relations to the mathematics of systems of integrable differential 
equations [7] and moduli space of Riemann surfaces [8, 9], make this a fascinating and novel 
(ab)use of quantum field theory. While condensed matter problems involve the thermal 
fluctuations of observable surfaces in 3D Euclidean space, for applications to high-energy 
physics the surfaces are worldsheets which enter indirectly through the weak coupling ex- 
pansions of relativistic string theories purporting to describe quantum gravity and confining 
gauge theories. Matrix models of relativistic strings should be defined in Minkowski space, 
their spectrum, or a self-consistent truncation of it, comprising string excitations. While 
knowledge of the fractal geometry of perturbative diagrams is sometimes useful, numerical 
computation of critical random surface properties in Euclidean space is no substitute for di- 
rect calculation of the relativistic string observables; for example, the almost uninvestigated 
issues of strongly coupled string theories would be otherwise neglected. 

I.Klebanov and the author have suggested [10] that discretised fight-cone quantisation 
(DLCQ), introduced with a view to calculating the bound-state spectrum of gauge theo- 
ries directly [11], may be very appropriately applied to matrix models also. This method, 
although primarily a numerical one, paints a clear physical picture of otherwise complex dy- 
namics. In the case of string theories, this means the size, shape, and energy of strings. The 
unfortunate circumstance that the simplest bosonic string theory, to which one would add 
more structure for a physically realistic object, is (expected to be) tachyonic poses no partic- 
ular problem to DLCQ of matrix models. On the contrary, one may analyse this pathology 
with some rigour and, as a result, set about curing it. 

The following section contains an elementary review of matrix models, light-cone quan- 
tisation for a c = 2 model being carried out in section 3. Section 4 describes the critical 
behaviour of the discretised version, while future avenues of development are discussed in 
section 5. 
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2. Matrix Models of Random Surfaces 



Consider an NxN hermitian matrix field (pabi^) in c dimensions subject to the following 
(Euclidean) action 



N 



The Feynman rules are those of ordinary 0^ field theory except that worldlines are double 
lines, each line carrying a "colour" index a. Colour is conserved along the propagator and 
at vertices on account of the global U(N) symmetry, — > Q^(f)Q, of the action; 

< <l)ab{x)(/)cd{y) >= ^adhc I TTT^^- > (2) 

< (t>ah{x)(t>cd{x)(t>ef{x) >= hc^de^faN ■ (3) 

The 1/N expansion of the theory described by (1) is a topological expansion in surfaces 
[1], the Feynman diagrams being understood as fishnet approximations drawn on continuous 
2-dimensional surfaces embedded in R^; 



J D0e-^^ -^J^^'' (;^) ^ J embeddings 



(4) 



V is the number of vertices, x the Euler number, and one integrates over all embeddings of 
the graph with rules (2) (3) . The discretised surface picture is formalised by considering 
the dual graphs [3, 4] that join centres of neighbouring loops, which for a 0^ theory specifies 
a simplicial triangulation. Each triangle carries unit intrinsic area and its centre coincides 
with the Feynman vertex; in particular this leaves the angle between neighbouring triangles 
unspecified in this model. The idea is then to tune the coupling A to a critical value Ac at 
which surfaces of large intrinsic area v are favourable. At this point one may be able to take 
a continuum limit for surfaces* 

★ There may be no scaling of the intrinsic geometry (curvatures, etc) even though the area scales - the 
c < 1 models arc examples of this - while at c > 1 there is even a danger, commonly attributed to 
tachyons, of non-scaling area. 
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The obscrvablcs are given by the Green's functions of the matrix field theory. A closed 
string state is a hole cut into the surfaces at some fixed time, to say. Such states are given 
by the singlet operators in the matrix model; 



Tr[<j>{xi)...<j>{xB)]t=to 



(5) 



is a B-hit string, the bits being embedded at xi, . . . ,xb respectively and forming a closed 
chain dual to the external legs of the Feynman diagrams. Also present are direct products 
of (5) (multi-string states) and non-singlets; the latter do not have an interpretation as 
closed strings and should be eliminated, either by performing a self-consistent truncation or 
effecting a dynamical decouphng (by gauging the U{N) for example [12]). The critical point 
Ac of large surfaces should manifest itself as singularities of Green's functions, such as the 
string propagator ~< Tr[0 • • • 0]Tr[0 • • • 0] >, and thus in particular the spectrum of string 
excitations (5) should exhibit some sort of critical behaviour. Applying DLGQ to the matrix 
field theory in Minkowski space, one derives the spectrum as a function of A and can search 
for such behaviour. 

The simplest non-trivial example to consider is c = 2, for which (1) models, in its 1/N 
expansion, a discrete version of the worldsheet action for the c — 2 non-critical bosonic string 



A, u, and T are renormalised parameters associated with A, and fi respectively. Stretch- 
ing energy of the worldsheet is governed by the propagator (2) , which specifies the probability 
amplitude for the separation of centres of neighbouring triangles. This exponential fall-off 
gives the Gaussian term in A plus non-renormalisable higher derivative terms. Naively the 
latter are irrelevant, but such reasoning assumes scaling of the worldsheet in some sense. The 
c — 2 theory (1) is rendered perturbatively finite by normal ordering, this being necessary 
in any case to eliminate graphs dual to pathological triangulations where two or more sides 
of the same triangle are identified. According to general arguments [14], the sum of graphs 
at a given order in 1/A^ in a UV finite theory is also finite for sufficiently small coupling 
constant; the 1/A^ expansion itself is only asymptotic. As A — > Ac one then approaches the 



[13], 




(6) 
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edge of the domain of convergence. The phenomenon is similar to, but certainly different 
from, the crossover to non-borel summability in the non-matrix field theory [N = 1). 



3. Light-Cone Quantisation 

Let us now rotate —>■ ix^ and find the relativistic string spectrum by light-cone 
quantisation. Defining light-cone variables — (a;^±a;^)/-\/2, with x~^ — x-, the light-cone 
energy and momentum — J dx~T~^^ are 



p-(x+) = I dx-Txi-iKt? ^03) . 



(7) 



The light-cone Hamiltonian P propagates a field configuration from one x"*" slice to another. 

Choosing a free field representation at x+ = say*. 



oo 





and imposing the canonical commmutation relations at equal x'^ lines, 

[(pijix'), d-(j)ki{x~)] = -S{x~ - x~)SiiSjk , (9) 
the modes a^- satisfy standard commutators; 

[aij{k+), al{k+)] = S{k+ - k+)6iiSjk ■ (10) 



energy quanta^ In the Fock space constructed using (10) , the orthonormal single closed- 



An important feature is that longitudinal momentum k'^ is positive semi-definite for positive 
energy quanta^ In the Fock 
string states are the singlets 

N-^/'^Tr[a\k+)---a\k+)]\0> , = ■ (H) 

If we let ^ oo the multi-string states can be neglected, since 1/A^ is the string coupling 
constant, and P^ propagates without splitting or joining strings. Thus we solve for the free 



B 



★ The symbol f is always imdcrstood to have purely quantum meaning and never acts on indices, 
t This positivity constraint results in the extra 1/2 in (9) , ensuring for example that [P^ , (t>] = d-(j). 
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string spectrum. (Non-singlcts states arc discarded by hand.) P+ and P can be simulta- 
neously diagonalised and so eigenfunctions of P~ with given will be some superposition 



P+ 



* = 5Z [ dki---dkBS{P+ -J2^i^-fB{ki,...,kB)N-^/^TT[a\ki)---a\kB)\0> . (12) 

B 



Explicitly one finds 



oo 







X 



J I ^ ^ ^ = {^jif't + k+)aik{k+)akj{k+) + al{k+)al.{k+)aij{k+ + k+)] 

(13) 

where repeated indices are summed over. For the free string theory each of the terms in 
(13) has a simple local action on the string. The mass term acts as a tensional energy 
~ fi l/ki and does not change the state. The cubic term can coalesce two neighbouring 
bits in the trace, or perform the inverse process, and changes both the length and momentum 
distribution (structure function) of the string. Since c = 2, there arc no transverse oscillations 
in the target space. Unlike c < 1 however, we are dealing with truely stringy degrees 
of freedom - an infinite number of particle fields - which supplement the centre of mass 
motion. 

Free string states will satisfy a relativistic dispersion relation 2P'^P~ — and by 
diagonalising P~ in the basis of states of a fixed P~^ we can find the spectrum of masses 
This diagonalisation cannot be performed analytically for the full Hamiltonian in general, so 
we must introduce a cutoff [11, 15], rendering the number of states of momentum P^ finite, 
and compute numerically. It is therefore important that we can neglect zero momentum 
modes a^(0) since we could in principle include arbitrary numbers of them without changing 
P"*". Such modes have infinite energy for 7^ 0, according to (13) ; however, they can be 
infinite in number and so a more careful analysis is required to lift this ambiguity. Indeed, it 
is believed by some that a proper constrained quantisation of these zero modes is necessary 
to describe spontaneous symmetry breaking and other non-perturbative effects in light-cone 
formalism. In those cases it is argued that the true vacuum, if one exists at all, is not the 
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Fock one, : P~ : |0 >= 0, but a condensate of a^(0) modes. But one should recall that we are 
interested in (1) only insofar as it generates random surfaces through the 1/N expansion and 
perturbation theory. Therefore we should always take the = vacuum, for which there is no 
condensate of zero modes, since this is the one with respect to which the planar diagrams are 
defined. For the c — 2 model, at each order in 1/N we work with a convergent perturbation 
expansion, the non-perturbative effects presumably manifesting themselves through the e~-^ 
corrections to the asymptotic expansion in 1/N. The latter are non-perturbative effects of 
string theory and specifying their details is equivalent to stating how one is going to stabilise 
the unbounded (f)^ theory (without disturbing the 1/N expansion), a question that will not 
be considered here. 

Another useful way of viewing these and other issues is to consider triangulations in 
light-cone perturbation theory. Indeed, for relativistic strings one could have set up the 
random surface expansion from this vantage point from the very beginning. Using Feynman's 
causality prescription on x'^ rather than the single-particle propagator becomes 

(oo \ 
r dp+ ^ ^+ _j(a;-p++x+/j/2p+) ^ f dP'^ g/ ^+^-i(x-p+-x+u/2p+) I ^ ^(^^) 
J 47rip+ J Amp+ j 2^11 

€ -OO / 

(14) 

The first two terms can be given the usual particle and anti-particle interpretation by viewing 
the negative energy (p~) states as propagating backwards in time (x"*"). The third term is 
a special contribution from p'^ — and corresponds in the dual diagram to the propagation 
of a zero momentum string-bit; neglecting this case, all particles and anti-particles move 
forward in x'^ carrying positive p'^. In early work [16], the quantisation procedure used was 
the one adopted here, where it was proved not only that ^'''-ordered and x°-ordered pertur- 
bation theory are equivalent, but also that the third term eventually does not contribute in 
non- vacuum diagrams. This seems to indicate once again that, provided we restrict to per- 
turbation theory of (1) , the zero modes can be ignored in computing the string propagator. 

4. Discretisation and Critical Behaviour 

The desired cut-off will be introduced by compactifying x~ and imposing periodic bound- 
ary conditions, Mij{x~) = Mij{x~ + L) [11]. Then the allowed momenta are labelled by 
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positive integers nf, 

kr = ?^,P- = ?^,f:n,=K. (15) 

1=1 

For fixed P^, removing the cut-off L ^ oo corresponds to sending K oo. The "harmonic 
resolution" K represents the total number of momentum units available to the string. The 
longest string has K bits of one unit*, the shortest one bit of K units, and in general the 
states can be labelled by the ordered partitions of K modulo cyclic permutations. Light-cone 
quantisation (8) - (13) may now be repeated for discrete variables A;''" — > n and one finds the 
mass relation 

op+p— \ 

= K(V -xT)- -—^ . (16) 

V is the discrete version of the mass term in (13) while T is the cubic term. For finite K 
the r.h.s. is a finite-dimensional symmetric matrix with real dimensionless elements - /i is 
the quantity of dimension mass^ which plays the role of string tension Ija' - which may 
be diagonalised as a function of the dimensionless parameter x] V is diagonal while T is 
off-diagonal. 

The following picture of the critical behaviour, supported by the numerical results [10, 
17], begins to emerge. In the light cone formalism of critical string theory the longitudinal 
momentum supported between two points on the string is proportional to the amount of 
cr-coordinate space between these points. We can adopt a similar co-ordinate system for the 
non-critical strings (5) . Indeed, fixing a particular bit as origin, we can define a positive 
scalar field on this cr-space hy X — Afe/Ao", where h is the distance, measured in number of 
bits, from the origin. For example, the zero mode J Xda is the intrinsic length of the string. 
As we remove the cutoff on the longitudinal momentum allowed for bits {K — > oo), and 
hence on discreteness of cr-space, the scalar field will generically take constant values almost 
everywhere in cr-space. We would like to be able to tune the theory to a critical point where 
the scalar field is in a long wavelength regime. In this case it would be somewhat similar to 
the Liouville mode of Poyakov's string [13]. What would this long wavelength regime mean 
for the spectrum? Firstly we would expect long string dominance; the expectation value 
of length for low-lying eigenstates would typically diverge (and we have been 



★ This sector alone represents what one might call "critical string theory" [24] 
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assuming that this is not distinct from Ac discussed earUer). Since, roughly speaking, each 
string-bit carries finite energy, we expect that |M^| ^ oo as a result. — * oo at x = 
for infinitely long strings, but as x ^ Xc, if only for consistency, we must see ^ —oo 
for the low-lying eigenstates if they are long. Indeed this will tend to happen to the lowest 
eigenvalues of any real symmetric matrix as one increases its dimension for sufficiently large 
off-diagonal elements. Only if the Hamiltonian is an explicitly bounded operator combination 
{H ~ O^O + const) can this instabihty be avoided in general. We might also expect to see 
a continuous spectrum at x — Xc if we compare with the Liouville theory results [18] 

2P+p-a' =pl + 4r-^, (17) 

but it has been difficult to confirm this numerically. Moreover the groundstate has finite 
negative mass squared in (17) , while the matrix model's is infinitely negative at x — Xc- 
It has been suggested [17] that perhaps fi should be renormalised to zero as a result, but a 
derivation of this requirement is still lacking. In any case we have the first direct demon- 
stration that the non-critical bosonic string with area action is tachyonic above c — 1. The 
use of an unbounded 0^ matrix potential does not a priori spell tachyons at any order in 
1/N in the expansion about = 0, but the string theory described at the critical coupling 
is nevertheless tachyonic.^ 

5. Future Directions 

In order to identify the critical behaviour at c = 2 more precisely it is useful to investigate 
the effects of adding an explicit polymerisation term [19] 

S = Se + J d^x ^(Tr[02])2 . (18) 

For sufficiently large g this worldsheet contact interaction seems to favour short strings. 
This is quite unlike the critical point at = and, assuming that there is a phase transition 
somewhere in between, casts doubt on branched polymer behaviour of the c — 2 matrix 

f The use of an unbounded potential is rather a symptom of the divergence of the 1 /N string perturbation 
expansion, as commented earUer. 
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model in Minkowski space; recall that this behaviour was identified in c > 1 Euclidean 
dimensions from numerical simulations [20] and combinatorial estimates [21] of dynamical 
triangulations at the critical point. Moreover simple polymerisation is not the only possibility 
in the Euclidean game. The phase diagram needs to be investigated in more detail before a 
clearer picture can be gained. 

The tachyon we found is expected to persist at c > 2. To regulate the transverse 
dimensions we can use a transverse lattice. To ehminate the zero modes associated with 
x±, which would otherwise obscure the stringy part of the spectrum, we can perform the 
Eguchi-Kawai (EK) compactification [22] to single links in each direction. The resulting 
field theories are much the same as the c = 2 one; namely, we deal with UV finite two- 
dimensional field theories with convergent perturbation expansions at given order in 1/N. 
If we use Hermitian matrix models (1) the EK reduction induces more general interactions 
Vef[{(f)) in the effective c — 2 potential [23]. Unfortunately V^ff(0) is not known explicitly and 
can only be calculated as an expansion in powers of the inverse transverse lattice spacing 
1/a. Truncating the expansion arbitrarily at some order, the resulting model will exhibit 
polymeric behaviour for sufficiently small a since the leading term is a contact interaction 
similar to that appended in (18) , the induced g being ~ 1/a^. At sufficiently large a the 
(unreduced) transverse lattice sites are obviously uncoupled and we should recover copies of 
the c = 2 model at each site. 

In order to deal with an exact EK reduced system, one can employ complex matrix 
models similar to the old Weingarten model [2] . Complex matrices M five on the links of a c- 
dimensional hypercubic lattice. For non-critical string theory we must use an action given by 
traces around all (oriented) loops of length 4 on the lattice [6] , which comprises the standard 
plaquette action plus zero area loops. Expansion of this theory in the coupling constant and 
1/A^ reproduces the random surfaces of a dynamical quadrangulation in c dimensions, the 
the probability distribution for neighbouring vertices {x, y} of the quadrangulation in this 
case being 



for orthonormal vectors p,. For c = 1 it has been proven that this gives the same answers 
at the critical point as using Feynman propagators [6]. To perform DLCQ we must take 
the naive continuum limit for two of the dimensions, which seems to produce an intractable 




(19) 



10 



two-dimensional kinetic term unfortunately. Instead we could use the Feynman propagator 
for these two continuous dimensions and study actions like [24, 25] 

J (fxTtidaMd'^M^ + fxMM^ + XiMMM^M^ + X2MM^MM^] . (20) 

This is a c = 3 model with EK reduction of the transverse direction - a single complex 
matrix M lies on the periodic link - particularly simple since it has no standard plaquette 
term. One may now study the DLCQ as a function of a. At sufficiently small a there should 
be a roughening transition from the c = 2 to the c — 3 phase. It should be noted however 
that the computational accuracy diminishes significantly for c > 2 due to the increase in 
degrees of freedom, each string bit now carrying at least one more Z2 variable - e.g. the real 
and imaginary parts of M. 

Even if reliable data on models such as (20) can be collected, we still expect tachyons. 
We must try other possibilities to eliminate them, such as introducing some dependence 
on extrinsic geometry, which is also of interest in condensed matter problems. This is 
certainly possible for the complex matrix models, at least on the lattice target space, since 
the orientation of simplices is rather manifest. At a more fundamental level, we should look 
for tachyon-free matrix models with long string dominance if we wish to describe interacting 
continuum strings. As indicated earlier, such would be a delicate theory; using a positive 
definite Hamiltonian tends to counteract precisely the requirement of the critical point - 
crossover to divergence of perturbation theory. Also the obvious fact should be stressed 
that, contrary to the popular trends of research, such a theory need not possess an attendant 
continuum worldsheet expansion to be a sensible string theory. 

For other applications of light-cone matrix models, such as to confining gauge theories, 
continuity of strings may not be so important. For example the interesting Regge-like tra- 
jectories found in two-dimensional gauged matrix models [12] are most probably the result 
of restriction to sectors of fixed (discrete) string length*. According to the suggestion made 
earlier, this freezes the Liouville-like zero mode and exposes the quasi-harmonic energy levels. 
While it is unrealistic to expect solution of the gauged models in higher dimensions, since 
they are more complicated than pure large- A?" QCD, the two-dimensional models may help 

* In rcf.[12] low- lying mass eigenstates tended to consist of strings of some given length for mysterious 
dynamical reasons. 
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us to understand more clearly the relationship between gluonic and fundamental strings. 
They describe a limit of higher dimensional gauge theory in which all transverse directions 
are very compact; xj^-independent transverse potentials Aj_{x~^ ,x~) play the role of matter 
(f) in the gauged c — 2 matrix model. 

Clearly there are many interesting questions in applications of string theory which may 
be addressed by the light-cone matrix models through analytic and numerical techniques. 
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